Ecology and Evolution 

Unique inhibitory cascade pattern of molars in canids 
contributing to their potential to evolutionary plasticity 
of diet 



Developmental origins that guide the evolution of dental morphology and den- 
tal formulae are fundamental subjects in mammalian evolution. In a previous 
study, a developmental model termed the inhibitory cascade model was estab- 
lished. This model could explain variations in relative molar sizes and loss of 
the lower third molars, which sometimes reflect diet, in murine rodents and 
other mammals. Here, I investigated the pattern of relative molar sizes (inhibi- 
tory cascade pattern) in canids, a taxon exhibiting a wide range of dietary hab- 
its. I found that interspecific variation in canid molars suggests a unique 
inhibitory cascade pattern that differs from that in murine rodents and other 
previously reported mammals, and that this variation reflects dietary habits. 
This unique variability in molars was also observed in individual variation in 
canid species. According to these observations, canid species have greater vari- 
ability in the relative sizes of first molars (carnassials), which are functionally 
important for dietary adaptation in the Carnivora. In conclusion, an inhibitory 
cascade that differs from that in murine rodents and other mammals may have 
contributed to diverse dietary patterns and to their parallel evolution in canids. 



Masakazu Asahara 

Graduate School of Sciences, Kyoto University, Kyoto, 606-8502, Japan 

Keywords Abstract 

Canidae, Carnivora, dental formulae, dental 
morphology, evolvability, inhibitory cascade. 

Correspondence 

Masakazu Asahara, Graduate School of 
Sciences, Kyoto University, Kyoto 606-8502, 
Japan. Tel:-i-81 75 753 7731; Fax:-i-81 75 753 
3276; E-mail: kamonomana@zoo.zool. 
kyoto-u.ac.jp 

Funding Information 

Funded by Mishima Kaium Memorial 
Foundation, Japan Society for the Promotion 
of Science, and Ministry of Education, 
Culture,Sports,Science and Technology in 
Japan 

Received: 11 October 2012; Revised: 25 
October 2012; Accepted: 1 November 2012 

Ecology and Evolution 2013; 3(2): 278-285 

doi: 10.1002/ece3.436 



Introduction 

Dental morphology and dental formulae are important 
taxonomic traits in mammals (Ungar 2010), and are also 
used for paleoecological and ecomorphological studies in 
mammals because these traits reflect dietary adaptations 
(Popowics 2003; Benton 2004; Friscia et al. 2007; Van 
Valkenburgh 2007). Patterns of adaptation are guided by 
the variability and evolvability of these traits (Klingenberg 
2005; Barton et al. 2007). In fact, evolvability, that is, sys- 
tems having variability, generating new variation (Wagner 
and Altenberg 1996; Kirschner and Gerhart 1998), is 
fundamental to the evolution of traits (Klingenberg 2005; 
Futuyma 2010). Therefore, developmental mechanisms 
that guide and constrain patterns of adaptation in dental 
morphology and dental formulae are crucial subjects for 
elucidating their proximate and ultimate factors, and the 
interactions between these factors during dental evolution 
(Kavanagh et al. 2007; Laland et al. 2011; Wilson 2011). 
Therefore, many recent studies have focused on the devel- 



opmental mechanism of the evolution of dental morphol- 
ogy and dental formulae, which guide and constrain the 
reflection of dietary adaptations (e.g. Kavanagh et al. 
2007; Polly 2007; Laffont et al. 2009; Renvoise et al. 2009; 
Harjunmaa et al. 2012; Wilson et al. 2012). A recent 
developmental study established a developmental model 
that can explain evolution of the relative sizes of lower 
molars in murine rodents (Kavanagh et al. 2007). This 
model, termed the inhibitory cascade model, explains the 
relative sizes of the lower molars (first, second, and third 
molars; Mj, M2, and M3, respectively) by the balance of 
inhibitor molecules from Mj tooth germ and activator 
molecules from mesenchyme during dental development. 
Inhibitor molecules inhibit the development of distal 
molars, whereas activator molecules activate their devel- 
opment. For example, greater inhibition generates a larger 
Ml and smaller M3 (Mj » M2 » M3), while lower inhi- 
bition and greater activation generate equal-sized molars 
(Ml = M2 = M3), and moderate inhibition and activation 
generate an intermediate condition (Mi > M2 > M3) 



278 



© 2012 The Author. Published by Blackwell Publishing Ltd. This is an open access article under the terms of the Creative 
Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original 
work is properly cited. 



M. Asahara 

(Kavanagh et al. 2007). This model can explain dental vari- 
ations that have resulted from dietary adaptations in mur- 
ine rodents (Kavanagh et al. 2007). Faunivorous murine 
species exhibit Mi » M2, and have lost M3. Conversely, 
herbivorous murine species have approximately equal- 
sized molars. For most mammals, from marsupials to vari- 
ous placental orders, the relative sizes of Mi, M2, and M3 
change sequentially and thus were explained by the model 
(Polly 2007). Several authors have investigated relative 
molar sizes in several mammalian taxa and have reported 
differences in the inhibitory cascade pattern between a 
number of taxa including murine rodents (Renvoise' et al. 
2009; Wilson et al. 2012). It has been noted that the vari- 
ability of a trait initiates the evolvability of that trait (Klin- 
genberg 2005; Barton et al. 2007; Wilson 2011). It is 
possible that the unique patterns of inhibitory cascade that 
guide variability in a particular taxon could promote the 
evolvability of typical molar patterns and, consequently, 
the evolvability of diet in that taxon. 

The order Carnivora, and particularly the family Cani- 
dae (canids), is one of the most diverse mammalian taxa 
in terms of dietary pattern, and exhibits parallel evolution 
in diet (Van Valkenburgh and Koepfli 1993; Friscia et al. 
2007; Goswami 2010; Sillero-Zubiri 2010). Similar to fau- 
nivorous murine rodents, several canid species have lost 
M3 (SiUero-Zubiri 2010); this loss is thought to be related 
to the enlargement of carnassial teeth (Mi) - a carnivo- 
rous adaptation for shearing flesh - and to the degenera- 
tion of molars M2 and M3 (HoUiday 2010). However, 
patterns of relative molar sizes and dietary adaptation in 
canids, and the relationship of these parameters to the 
inhibitory cascade model are still not clear. In this study, 
my primary objective was to elucidate patterns of inter- 
specific variation in the relative sizes of lower molars 
in canids, and to determine the relationship of this varia- 
tion to the inhibitory cascade model and to dietary 
adaptations. 

A second objective was to elucidate the variability within 
the species in canids that guides evolutionary patterns 
(Klingenberg 2005). To achieve these objectives, I investi- 
gated individual variation in relative molar sizes as an indi- 
cation of variability in this parameter at the intraspecific 
level (Klingenberg 2005). In addition, I investigated indi- 
vidual variation in the number of teeth, as oligodonty 
(missing teeth) is considered a transitional stage in the evo- 
lution of dental formulae (e.g. Ohtaishi 1986; Giannini and 
Simmons 2007). An earlier experimental study examining 
mouse development found that the number of molars was 
affected by the inhibitory cascade (Kavanagh et al. 2007). 
Therefore, I also compared individual variation in rela- 
tive molar sizes and number of molars (i.e. congenital 
missing of M3), to consider the inhibitory cascade and the 
evolutionary process of M3 loss in canids. 



Unique Inhibitory Cascade in Canid IVIolars 

Material and Methods 

I examined 320 specimens from 27 species of canids (Cani- 
dae, Carnivora, Mammalia) (Table 1). All species were 
examined to clarify evolutionary patterns in relative molar 
sizes. The dietary pattern of each species was categorized as 
carnivorous (primarily eating mammalian flesh), omnivo- 
rous (eating various foods, with neither mammalian flesh 
nor insects comprising >50% of the diet), or insectivorous 
(primarily eating insects) using information from the liter- 
ature (Sillero-Zubiri 2010). In order to estimate variability 
in molars, I examined individual variation in relative molar 
sizes. For this purpose, individual variation within seven 
species was examined whereby I measured >15 individuals 
from each species (Table 1, 2). In addition, I examined 
individual variation in the presence or absence of M3 in 
Vulpes lagopus and Nyctereutes procyonoides to clarify 
whether individual variation and missing of M3 are 
explained by the inhibitory cascade model. The specimens 
of N. procyonoides examined were those deposited in Kyoto 
University Museum, Kyoto University, Japan, which had 
been collected from a small island, Chiburi Island, Shimane 
Prefecture, Japan. Specimens of the other species were 
those deposited in the Department of Mammalogy, 
American Museum of Natural History, USA, which had 
been collected from large areas. I measured the size of each 
molar as the projected area in photos taken from the occlu- 
sal view using ImageJ software (NIH, Bethesda, MD), and 
compared the relative molar sizes in the morphospace: M2 
size/Ml size versus M3 size/Mi size (abbreviated as M2/M1 
vs. M3/M1) (Kavanagh et al. 2007). Any given point in 
morphospace represents the relative sizes of the three 
molars of a particular species or individual. I plotted the 
average values for each species to describe interspecific vari- 
ation, and plotted each individual to describe individual 
variation. Reduced major axis (RMA) regressions were per- 
formed on these plots after performing Anderson-Darling 
normality test. I used M2/M1 as an index of activation ver- 
sus inhibition during molar development. M2/M 1 scores 
between carnivorous and omnivorous species were com- 
pared using the Mann-Whitney U test. Further, for V. lag- 
opus and M procyonoides, M2/M1 scores were compared 
between normal individuals and individuals that were miss- 
ing M3 on one or both sides. Statistical analyses were per- 
formed using Minitab 14 (Minitab, Inc., PA), and RMA 
regressions were performed using PAST (Hammer et al. 
2001). Several studies have utilized multiple regressions to 
elucidate how absolute molar sizes affect one another (Ren- 
voise' et al. 2009; WUson et al. 2012). However, this 
method tends to reflect variability in the absolute size of 
Ml, and activation versus inhibition patterns can become 
obscured. Therefore, I focused on relative molar sizes, that 
is, the inhibitory cascade. 
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Table 1. Species examined in this study, and their molar ratios. 
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Number 


Species 


Diet 


W 


M2/M1 ± SD 


I\/!3/M1 ± SD 


1 


Atelocynus microtis 


Omnivorous 


3 


0.54 ± 0.04 


0.15 ± 0.04 


2 


Canis adustus 


Omnivorous 


4 


0.51 ± 0.12 


0.17 ± 0.05 


3 


Canis aureus 


Omnivorous 


11 


0.42 ± 0.04 


0.1 1 ± 0.02 


4 


Canis latrans 


Carnivorous 


51 


0.37 ± 0.03 


0.09 ± 0.02 


5 


Canis lupus 


Carnivorous 


28 


0.31 ± 0.03 


0.09 ± 0.01 


6 


Canis mesomelas 


Omnivorous 


20 


0.36 ± 0.02 


0.10 ± 0.02 


7 


Cerdocyon thous 


Omnivorous 


5 


0.54 ± 0.06 


0.16 ± 0.03 


8 


Chrysocyon brachyurus 


Omnivorous 


3 


0.46 ± 0.01 


0.18 ± 0.01 


9 


Cuon aipinus 


Carnivorous 


2 


0.25 ± 0.00 


0.00 ± 0.00 


10 


Lycalopex cuipaeus 


Omnivorous 


10 


0.41 ± 0.03 


0.1 1 ± 0.02 


11 


Lycalopex griseus 


Omnivorous 


15 


0.50 ± 0.05 


0.13 ± 0.02 


12 


Lycalopex gymnocercus 


Omnivorous 


9 


0.53 ± 0.02 


0.13 ± 0.01 


13 


Lycalopex sechurae 


Omnivorous 


5 


0.53 ± 0.04 


0.16 ± 0.01 


14 


Lycalopex vetulus 


Insectivorous 


6 


0.71 ± 0.12 


0.25 ± 0.1 1 


15 


Lycaon pictus 


Carnivorous 


7 


0.32 ± 0.02 


0.07 ± 0.02 


16 


Nyctereutes procyonoides 


Omnivorous 


44 


0.47 ± 0.03 


0.07 ± 0.04 


17 


Otocyon megatotis 


Insectivorous 


7 


0.97 ± 0.05 


0.82 ± 0.06 


18 


Speothos venaticus 


Carnivorous 


4 


0.17 ± 0.04 


0.00 ± 0.00 


19 


Urocyon cinereoargenteus 


Omnivorous 


31 


0.53 ± 0.03 


0.14 ± 0.03 


20 


Vulpes bengalensis 


Omnivorous 


2 


0.61 ± 0.03 


0.23 ± 0.04 


21 


Vulpes chama 


Omnivorous 


1 


0.62 ± 0.00 


0.20 ± 0.00 


22 


Vulpes lagopus 


Carnivorous 


31 


0.32 ± 0.04 


0.07 ± 0.04 


23 


Vulpes macrotis 


Omnivorous 


4 


0.39 ± 0.05 


0.10 ± 0.00 


24 


Vulpes pallida 


Omnivorous 


1 


0.72 ± 0.00 


0.23 ± 0.00 


25 


Vulpes velox 


Omnivorous 


7 


0.39 ± 0.01 


0.09 ± 0.01 


26 


Vulpes vulpes 


Carnivorous 


3 


0.35 ± 0.01 


0.10 ± 0.00 


27 


Vulpes zerda 


Omnivorous 


6 


0.58 ± 0.04 


0.16 ± 0.03 



Results 

As a result of interspecific variation, plots of the molar 
ratios of all species in morphospace indicated that relative 
molar sizes changed sequentially (i.e. Mi > M2 > M3; 
plots are in the white zone in Fig. 1). These results are in 
agreement with the consensus area of the inhibitory 
cascade model suggested by PoUy (2007). Interspecific 
variation in relative molar sizes among the majority of 
canid species, excluding Otocyon megalotis, exhibited a 
pattern that differed in slope from the variation observed 
in murine rodents (Kavanagh et al. 2007). That is, the 
area in which O. megalotis was plotted indicated that it 
had similar sized molars (Fig. 1). Including O. megalotis 
data in M3/M1 scores caused the assumption of normality 
to be violated; therefore, this species was excluded from 
interspecific regression analysis. The pattern of interspe- 
cific variation revealed a correlation between M2/M1 and 
M3/M1 (P < 0.001; Table 2). The slope of the variation 
among murine rodents was higher than, and outside of 
the 95% confidence interval of, that of canids (Table 2), 
indicating that the cascade patterns of canids and murine 
rodents are significantly different (Fig. 1, Table 2). 
According to this pattern, the relative sizes of Mi in 



canids varied greatly in relation to that of murine 
rodents. Consequently, the maximum relative size of Mi 
occupied 80% of the total molar row in canids, but only 
66% of that in murine rodents (Fig. 1). Carnivorous spe- 
cies tended to have lower M2/M1 scores than omnivorous 
species (two-sided U test, W= 21.0, P < 0.001). That is, 
some carnivorous species have much larger Mi and smal- 
ler M2 and M3, than omnivorous species. Notably, species 
at the extremes of the distribution, that is, Cuon aipinus 
and Speothos venaticus, have lost M3. These differences in 
relative molar sizes and loss of M3 have evolved in paral- 
lel within many clades of canids (Fig. 3). 

As the result of individual variation, normality of the 
M3/M1 data was not observed for V. lagopus and 
N. procyonoides, species in which individuals were missing 
M3. When individuals with dental anomalies were 
excluded, the data for all species were normally distrib- 
uted (P < 0.05), and regression analyses were performed. 
In the M2/M1 versus M3/M1 morphospace, individual 
variations were correlated in most of the species 
(P < 0.05; Table 2), with the exception of Canis mesom- 
elas {P = 0.06; Table 2). That is, individuals with rela- 
tively larger Mj tended to have relatively smaller M2 and 
M3, and vice versa (Fig. 2). Ten individuals of V. lagopus 



280 



© 2012 The Author. Published by Blackwel! Publishing Ltd. 



M. Asahara 



Unique Inhibitory Cascade in Canid IVIolars 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
M2/M1 



and 25 individuals of N. procyonoides were missing M3 
on one or both sides (32% and 56%, respectively). There 
was no evidence of concrescence, and all cases of missing 
teeth were considered to be congenital. In both V. lag- 
opus and N. procyonoides, individuals missing one or two 
M3 had lower scores for M2/M1 (i.e. greater inhibition) 
than normal individuals (one-sided U test, W = 169.0 
and 509.0, respectively, P < 0.05) (Fig. 2). Individuals 
with relatively larger Mj and smaller M2 tended to be 
missing M3. 

Discussion 




80% 

Carnivorous Omnivorous Insectivorous 



Figure 1. Variation in relative molar sizes among canid species. 
Interspecific variation in canids, excluding Otocyon megalotis (red 
line), differed from that in murine rodents (blue line). Colors and 
shapes indicate diet of a given species (red square: carnivorous, blue 
circle: omnivorous, green triangle: insectivorous) (Table 1). (a) Species 
plots with standard deviations (SD). Numbers indicate species as in 
Table 1. (b) Dietary patterns with SD. Occlusal view of molar rows of 
species for each diet in canids and murine rodents. Illustration of 
murine rodents after Kavanagh et al. (2007). Maximum relative M, 
size reaches 80% of the total molar row in canids, and 66% in 
murine rodents. 



The inhibitory cascade pattern in canids and 
its relationship to diet 

Sequential changes in relative molar sizes (i.e. Mj > M2 > M3; 
plots are in the white zone in Fig. 1 ) and correlation between 
M2/M1 and M3/M1 have been considered as evidence that 
relative molar sizes are regulated by an inhibitory cascade, 
indicating that there are single mechanisms that inhibit 
distal molars (Kavanagh et al. 2007; Polly 2007). Therefore, 
variation in relative molar sizes in canid species is also regu- 
lated by an inhibitory cascade. In this study, carnivorous 
species tended to have lower M2/M1 scores (or relatively 
larger Mi) than omnivorous species. Moreover, the two 
carnivorous species having the smallest M2 in relation to Mj 
have lost M3. Carnivorous species exhibited the pattern 
Ml » M2 » M3, but omnivorous species exhibited 
Ml > M2 > M3. As the number of insectivorous species was 
limited, they could not be analyzed statistically; however, 
these species tended to have more equal-sized molars. Thus, 
the inhibitory cascade reflects dietary adaptation in canid 
molars (Fig. 1). The relationships among inhibitory cas- 
cade, relative molar sizes, and dietary adaptation in canids 
are similar to those in murine rodents (Kavanagh et al. 
2007). However, the patterns of adaptation differ between 
canids and murine rodents. For example, insectivorous 
canids and herbivorous murines have equal-sized molars. 



Table 2. Regression results (RMA) of the M2/M1 versus M3/M1 morphospace, showing confidence intervals (CI). 



Types of variation Slope CI max CI min Intercept CI max CI min r P N 



Inhibitory cascade model 




2.00 






-1.00 












Canidae (with 3 molars) 


Interspecific 


0.45 


0.515 


0.376 


-0.08 


-0.037 


-0.104 


0.91 


0.000 


24 


Canidae (without 0. megalotis) 


Interspecific 


0.48 


0.537 


0.412 


-0.09 


-0.057 


-0.119 


0.93 


0.000 


26 


Canidae (on diet) 


Diet 


0.48 


0.538 


0.438 


-0.09 


-0.061 


-0.111 


0.99 


0.035 


3 


Canis latrans 


Individual 


0.52 


0.632 


0.380 


-0.10 


-0.041 


-0.138 


0.41 


0.003 


51 


Canis lupus 


Individual 


0.49 


0.606 


0.267 


-0.07 


0.003 


-0.101 


0.69 


0.000 


28 


Canis mesomelas 


Individual 


0.71 


2.059 


0.240 


-0.16 


0.011 


-0.649 


0.43 


0.060 


20 


Lycalopex griseus 


Individual 


0.46 


0.594 


0.184 


-0.09 


0.036 


-0.163 


0.60 


0.018 


15 


Urocyon cinereoargenteus 


Individual 


0.94 


1.150 


0.676 


-0.36 


-0.220 


-0.466 


0.54 


0.002 


31 


Vulpes lagopus (with 3 molars) 


Individual 


0.73 


0.946 


0.425 


-0.16 


-0.052 


-0.233 


0.78 


0.001 


15 


Nyctereutes procyonoides (with 3 molars) 


Individual (small island) 


0.82 


2.096 


0.560 


-0.32 


-0.186 


-0.984 


0.55 


0.036 


26 
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(a) 

Normal 
Missing 



0.15 



0.10 



0.05 



0.00 



0.25 



0.30 0.35 
IVI2/M1 



0.40 



(b) 

Normal 
IVIissing . 




0.40 



0.45 0.50 
I\/I2/M1 



0.55 



Figure 2. Individual variation in relative molar sizes in (a) Vulpes lagopus and (b) Nyctereutes procyonoides, indicating presence or absence of M3 
by color (black: missing on both sides; dark gray: missing on one side; right gray: normal). Within species, M2/M1 scores differed significantly 
between normal individuals and individuals missing IVI3 (indicated by asterisks). 



and carnivorous canids and faunivorous murines (eating 
animals including insects) have relatively larger Mi. This 
may be due to the difference in absolute body size 
among the species. Insects are sufficiently large prey for 
murines, and these mammals need to concentrate their 
masticatory function on one major tooth. Similarly, mam- 
malian flesh is sufficiently large to be accommodated by 
canid molars. In contrast, insects are small food items for 
canids, and canids require a long molar row with equal- 
sized teeth in order to chew a number of insects at once 
(Ungar2010). 

Among the species examined, only O. megalotis had 
equal-sized molars, and was located distantly from the 
other canids in morphospace (Fig. 1). This may be related 
to the unique characteristic in O. megalotis of having four 
lower molars (Sillero-Zubiri 2010; Ungar 2010). Interspe- 
cific variation in the other canids exhibited a unique 
pattern of molar ratios that differs from that in murine 
rodents (Kavanagh et al. 2007), indicating a difference in 
the inhibitory cascade. Such differences have been 
reported in previous studies; however, the slope of the 
difference in canids was lower than that in any previously 
reported taxa (canids, interspecific: 0.48; canids, individ- 
ual: 0.46-0.94; murine and arvicoline rodents and South 
American ungulates: 1.17-2.15) (Kavanagh et al. 2007; 
Renvoise' et al. 2009; Wilson et al. 2012). In mouse exper- 
iments, all inhibition molecules were eliminated and 
interspecific variation in murine rodents was identical to 
observed variation in molar proportions (Kavanagh et al. 
2007). However, diffusion patterns may differ between 
inhibition molecules. It is possible that particular mole- 
cules with low diffusion efficiency have high evolvability 
and generate unique slopes in canids; however this is not 
yet clear. 



Variability in relative molar sizes and loss of 
M3 in canids 

The results of individual variation clearly showed the corre- 
lations between M2/M1 and M3/M1 (Table 2), and different 
M2/M1 scores between individuals having M3 vs. individu- 
als in which M3 was missing (Fig. 2), indicating that greater 
inhibition results in smaller distal molars and/or a loss of 
M3. Although the correlation coefficient r was not high 
between individuals (Table 2), individual variation reflects 
a large number of environmental factors. Therefore, a sig- 
nificant correlation indicates that individual variation 
reflects an inhibitory cascade. Individual variation indicates 
the variability within a species (Klingenberg 2005); there- 
fore, canid species differ from murine rodents in terms of 
variability, that is, they have lower slopes (Table 2). These 
patterns of variability are likely to be the source of the 
unique interspecific variation observed in canids. 

The results of interspecific and individual variation in 
M3 loss indicate that M3 loss in canids must be generated 
by greater inhibition during evolution. Individual plots 
for V. lagopus provided a good illustration of interspecific 
variation; V. lagopus individuals with missing M3 were 
plotted near C. alpinus, a species that has lost M3, 
whereas individuals with normal dentition were plotted 
near species that retain M3 (e.g. Vulpes vulpes, Vulpes 
macrotis) (Figs. 1 and 2). Therefore, V. lagopus may be in 
a so-called 'transitional stage' of evolution of dental 
formulae, reflecting greater inhibition of the inhibitory 
cascade along the trajectory of carnivorous adaptation 
(Fig. 3). Nevertheless, the N. procyonoides population 
exhibited a high frequency of missing M3 despite 
relatively lower inhibition in relation to V. lagopus or 
other canids that retain M3 (Figs. 1 and 2). This 
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I » I 23; Vulpes macrotis 

h*H 26: Vulpes vulpes 
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I — • — 1']9: Urocyon clnereoargenteus 

I • 1 1 : Atelocynus microtis 

' — • — I 7: Cerdocyon thous 

I • 1 1 1 ; Lycalopex griseus 

I % i12: Lycalopex gymnacercus 

■ ' — • — ' 1 3: Lycalopex sechurae 

*8; Chrysocyon bracchyrus 

18: Speothos venatlcus 

I • 1 2; Canis adustus 

I — • — I 3; Canis aureus 

- I ■ I 5; Canis lupus 

I ■ 1 4; Canis latrans 

- ■ 9: Cuon alplnus 

- - i-Bh 15: Lycaon plctus 



I 17; Otocyon megalotis 



0.1 0.2 0.3 0.4 



0.5 0.6 
M2/M1 



0.7 



0.9 1.0 



Figure 3. M2/M, scores for each species, and their phyiogenetic relationships from Bardeleben et al. (200B). Colors and figures indicate diet as in 
Fig. 1. Numbers indicate species as in Table 1. 



N. procyonoides population has probably been affected by 
inbreeding depression, as the population examined was 
introduced from the mainland to a small island (Saeki 
2009). The molar proportion of this population over- 
lapped with that of the mainland populations that retain 
M3 (personal observation). Geographical isolation and 
fixation of series of genes that relate to the inhibitory cas- 
cade or other mechanisms could also be an important 
process in the evolution of dental formulae (Asahara 
et al. 2012). The fact that no V. lagopus and N. procyono- 
ides individuals exhibited M3/M1 scores of lower than 
0.03 and 0.05, respectively, may relate to additional 
mechanisms for the regulation of M3 development and a 
possible threshold for M3 development or loss; teeth 
germ which are smaller than some threshold at the criti- 
cal stages cannot continue to develop into mature teeth 
(Gruneberg 1951; Wolsan 1989; Szuma 2003). 

Functional consequences of the unique 
inhibitory cascade pattern in canids 
contributing to the evolvability of diet 

The patterns of interspecific variation shown here are 
indicative of a unique inhibitory cascade pattern with less 
steep regression lines in the morphospace (Fig. 1, 
Table 2) than any other previously reported mammals 
(Kavanagh et al. 2007; Renvoise' et al. 2009; Wilson et al. 
2012). Guided by this pattern, the change in relative size 
of Ml has been amplified in canids (e.g. Mi comprises 
>80% of the total molar surface in S. venaticus, whereas 
the maximum proportion occupied by Mi in murine 
rodents is 66%; Fig. 1). My analysis is based on the two- 
dimensional occlusal surface, as used in previous studies 
(Kavanagh et al. 2007; Renvoise' et al. 2009; Wilson et al. 
2012). However, the canid Mi is a high cuspid tooth in 



relation to other canid molars, or all molars of rodents; 
therefore, if analysis is based on the three-dimensional 
tooth volume, the change in the canid Mi must become 
further amplified. 

Canids have evolved different proportions among func- 
tionally distinct parts of their lower molars, that is, the 
shearing surface (trigonid of Mi), which is important for 
a carnivorous diet, and the grinding surface (talonid of 
Ml, M2, and M3), which is important for omnivorous 
and insectivorous diets (Van Valkenburgh and Koepfli 
1993; Friscia et al. 2007). Therefore, the particular pattern 
of inhibitory cascade (with amplified change in MJ 
would contribute to dramatic changes in the proportion 
of shearing surface in Mi versus grinding surface in M2 
and M3 (i.e. dramatic changes in function) over the 
course of evolution. This dramatic change is regulated by 
a single developmental mechanism, the inhibitory cascade. 
Therefore, the molars of canids can readily evolve to 
adapt to a carnivorous, omnivorous, or insectivorous diet, 
and canids thus have the potential for evolutionary plas- 
ticity in their diet. These patterns of variation must have 
contributed to the diversity of dietary patterns and their 
parallel evolution among canids (Fig. 3) (Goswami 2010; 
SiUero-Zubiri 2010), and to the short- time divergence 
and diversity of dietary patterns in Lycalopex species (Pe- 
rini et al. 2009). PoUy (2007) inferred the existence of 
inhibitory cascade regulation across all mammals. In addi- 
tion, previous studies (Renvoise' et al. 2009; Wilson et al. 
2012), and my results, suggest that patterns of the inhibi- 
tory cascade can differ among taxa. Moreover, I suggest 
that these different inhibitory cascade patterns have con- 
tributed to different evolvability and diversity of diet 
among taxa. That is, clade-specific modification in develop- 
mental mechanisms could have promoted the capacity for 
dietary adaptation; that is, the dynamics of proximate and 
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ultimate factors. Investigation of inhibitory cascade patterns 
in other mammals will further our understEinding of these 
evolutionary dynamics. 
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